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Abstract 

We report studies of —>■ 7r+7r~7r°, using 78 fb~^ of data collected at the T(45) resonance 
with the Belle detector at the KEKB asymmetric e'^e^ storage ring. We measure the branching 
fraction for B^ ^ p^Tr^ to be B {B^ ^ p^n^) = (29.ll|S|(stat) ± 4.0(syst)) x IQ-^, and find 
an untagged charge asymmetry A = — 0.38lQ2i(stat)^Qg5(syst). We find the first evidence for 
B^ p^TT^ with 3.1(7 statistical significance and with a branching fraction of B (^B^ — > p^ir^^ = 
(6.0l||(stat) ± 1.2(syst)) x IQ-^. 

PACS numbers: 
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There is a large amount of interest in 'k~^'k~7t^ decays Ij. The three body final 

state is expected to be dominated by the quasi- two body decays P^^r^, the branching 

fraction of which has been previously measured to be B {B^ p^vr^) = (20.8j]Q3^3^) x 10~^ 
by Belle j2[, from a data sample of 29.4 fb~^, corresponding to 31.9 x 10^ BB pairs. Recently 
BaBar have announced a preliminary measurement of B (5° — > p^'^vr^) = (22.6 ± 1.8 ± 2.2) x 
10~^ [s^, from a data sample corresponding to 89 x 10^ BB pairs. Additionally, they have 
performed a time-dependent analysis on the quasi-two body signal candidates. Such an 
analysis can obtain information about the Unitary Triangle angle 02 which may be 
useful to help interpret the results of time-dependent B^ tt+tt" analyses 

In addition to p^'^tt^, the tt+tt^tt'^ final state can be accessed via B^ p^vr" decay. 
Analogous to B^ tt^tt" amongst the vrvr final states, B {B^ — > p°7r°) can be used to limit 
the possible contribution from penguin diagrams. Note that the penguin pollution in p^ir^ 
is expected to be smaller than that in tt+tt"; a hypothesis which is supported by comparing 
the ratio B (5° ^ p^K^) /B (5° ^ p^yr^) to i3 (5° n^K^) /B (5° ^ tt^tt^). Currently 
the best upper limit is B {B^ — > p^ir^) < 5.3 x 10~^ whilst most theoretical estimates of 
this branching fraction are O (10~^) or lower. 

An unambiguous measurement of 02 can, in principle, be made from the Dalitz plot of 
B^ Tr+Tr^Tr*^ [7J. It is also possible to observe direct CP violation from a population 
asymmetry in the untagged Dalitz plot distribution Furthermore, other resonant contri- 
butions to the 7r+7r~7r° final state are possible; recently there has been particular theoretical 
interest in 5° ^ an'^ Q. 

With very large statistics, and sophisticated analysis techniques, it would be possible to 
address all these open questions regarding tt+tt^tt^, using a time-dependent Dalitz 

plot analysis. While we cannot as yet achieve this ambitious goal, the results presented 
here represent milestones towards this objective. We present an updated measurement of 
the B^ — >■ p^tt"^ branching fraction, using event selection that can provide a sample for a 
quasi-two body time-dependent analysis. We also present the first evidence for B^ —>■ p°7r°. 

The analysis is based on a 78 fb~^ data sample containing 85 x 10^ B meson pairs 
collected with the Belle detector at the KEKB asymmetric-energy e'^e~ collider KEKB 
operates at the T(4S') resonance {^/s = 10.58 GeV) with a peak luminosity that exceeds 
1 X 10^'' cm~^s~^ 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of-fiight scin- 
tillation counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals 
(ECL) located inside a super- conducting solenoid coil that provides a 1.5 T magnetic field. 
An iron fiux-return located outside of the coil is instrumented to detect Ki mesons and to 
identify muons (KLM). The detector is described in detail elsewhere jll| . 

Charged tracks are required to originate from the interaction point and have trans- 
verse momenta greater than 100 MeV/c. To identify tracks as charged pions, we com- 
bine dE/dx information from the CDC, pulse height information from the ACC and timing 
information from the TOF into pion/kaon likelihood variables C{tt/K). We then require 
C{tt) / {C{tt) + C{K)) > 0.6. Additionally, tracks which are consistent with an electron 
hypothesis are rejected. 

Neutral pion candidates are reconstructed from photon pairs. Photon candidates are 
selected with a mininum energy requirement of 50 MeV in the barrel region of the ECL, 
and 100 MeV in its endcap. The vr" candidates are required to have momenta greater than 
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200 MeV in the lab frame, and a 77 invariant mass in the range 0.118 < M^^/GeV/c^ < 

0.150. Additionally, we require cos(^J[gj) < 0.95, where ^Jgj is defined as the angle between 

one photon's fiight direction in the vr" rest frame and the fiight direction of 7r° with respect 
to the lab frame, and make a loose requirement on the of a 7r° mass-constrained fit of 77. 

B candidates are selected using two kinematic variables: the beam-constrained mass 
Mf,c = \/ E^eam ~ ^^"^ energy difference AE = Eb — -Ebeam- Here, Eb and Pg 
are the reconstructed energy and momentum of the B candidate in the center of mass 
(CM) frame, and -Ebeam is the average beam energy in the CM frame. B candidates with 
Mbc > 5.2 GeV/c^ and -0.30 < AE/GeV < 0.20 are selected. We further define the signal 
regions: Mbc > 5.27 GeV/c^ and -0.10 < AE/GeV < 0.08. 

To select p^vr^ from 3-body tt+tt^tt^, we select candidates with an invariant tttc^ mass in 
the range \M^^o — Mp\ < 0.20 GeV/c^, and p helicity 9^^^, defined as the angle between the 
charged pion direction in the p rest frame and the p direction in the B rest frame l3|, in 
the range |cos6'(Jg[| > 0.5. 

The dominant background to B^ —>■ tt^ti^tt^ comes from continuum events, e"'"e" —>■ qq 
(g = u,d,s,c). Since these tend to be jet-like, whilst BB events tend to be spherical, 
we use event shape variables to discriminate between the two. We combine five modified 
Fox- Wolfram moments 0] into a Fisher discriminant; the coefficients are then tuned to 
maximize the separation between signal and continuum events. We further define 6'^ as 
the angle of the reconstructed B candidate with respect to the beam direction. Signal 
events have a distribution proportional to sin^(6'B), whilst continuum events are flatly dis- 
tributed in cos{6b)- We combine the output of the Fisher discriminant with cos(^b) into 
signal/background likelihood variables, Cs/b- We find the optimum selection requirement 
by maximizing S/ y/S + B, where S and B are respectively the expected numbers of signal 
and background events in the signal region. We use our measured branching fraction of 
B^ p'^Ti^ [5] as input, and find the optimum requirement is Cs/ {Cg + C-b) > 0.8. 

If more than one candidate remains in any event, that with the smallest Xvtx + X^o is 
selected, where Xvtx the of ^ vertex-constrained fit of vr+vr", and x^o is that from a 7r° 
mass-constrained fit of 77. 

We obtain the signal yield using a binned fit to the AE distribution, and cross-check the 
result by fitting the Mbc distribution. When fitting one variable, candidates are required 
to be in the signal region of the other. The signal probability density functions (PDFs) are 
obtained from Monte Carlo (MC); a Crystal Ball [15] lineshape plus a Gaussian is used for 
AE, whilst the Mbc distribution is described by a Gaussian. The Gaussian in the AE PDF 
accounts for poorly reconstructed low momentum neutral pions. The AE width is calibrated 
using an inclusive D* sample {D*~^ D^tt"*", — > K'tt'^tt^) whilst the AE and Mbc peak 
positions are adjusted according to a data sample of B^ D^tt~^ with — >■ K'n^ir^. 

The dominant background is from continuum events. The AE distribution for these 
events is described by a Chebyshev polynomial, whilst the Mbc shape is given by the AR- 
GUS function fl6'|. The parameters of these functions are determined from fitting a large 
continuum MC sample. 

Background is also possible from generic 6 — > c transitions; in this case the shape is hard 
to describe by a functional form and so we use a smoothed histogram. The AE distribution 
of background from the charmless decay B^ — >■ p'^p^ has a similar shape to the generic 
6 ^ c, so we combine these components, with the relative normalization fixed according to 
our recent measurement of B {B~^ — > p'^p^) allow the overall normalization to fioat 

in the AE fit. In the Mbc fit, these backgrounds cannot be distinguished from signal, and 
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FIG. 1: Results of fits to (left) IS.E and (right) Mbc- Contributions from signal (horizontally 
hatched), continuum (dashed line) and b ^ c + pp (diagonally hatched; denoted as BB in the 
legend) can be seen. The contributions from pK and other rare decays are small and hard to make 
out. The sums of contributions are shown as solid lines. 



the normalization is fixed. 

There are other possible backgrounds from charmless B decays, with distinctive AE 
shapes. p^K^ has a similar shape to the signal, but a shifted peak due to the 

misidentification of the kaon as a pion. The normalization of this component is fixed ac- 
cording to our recent measurement 18]. Contributions from , hh and hn*^ {h = tt^, K^) 
final states are scaled according to the most recent measurements of their branching frac- 
tions 0, 0, 2^, then combined into a smoothed histogram. The normalization of this 
component is then fixed in the fit. 

The results of the fits to AE and Mbc are shown in Fig. ^ From the AE fit we find a 
yield of 257.91 

43^2 signal events, with a significance of 6.3a. 

The significance is defined as ■>/— 2 ln(£o/'^max)5 where £max (-^o) denotes the likelihood 
with the signal yield at its nominal value (fixed to zero). 

From the Mbc fit we find a signal yield of 177.7l24o- Taking into account the different 
efficiencies of the AE/Mi^c signal region selections, this result is consistent with the AE 
yield. 

To check the events in the signal peak are B^ p^n^ events, and not from some other 
contribution to the tt+tt^tt^ final state, we relax the M^^t^o and cos^^^j criteria in turn, and 
perform the AE fit in bins of M^^o, and in bins of cos^(^gj. The resulting distributions are 
shown in Fig. |21 Clearly, the signal is consistent with being entirely due to B^ p^ir^ 
decays. 

To extract the branching fraction, we measure the reconstruction efficiency from MC 
and correct for a discrepancy between data and MC in the pion identification requirement. 
This correction is obtained from an inclusive D* control sample {D*~^ D^ir^, — > 
K~7^^), and is applied in bins of track momentum and polar angle. After this correction, 
the reconstruction efficiency is 10.4%. 

We calculate systematic errors from the following sources: PDF shapes (by varying 
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FIG. 2: Yield of signal events from IS.E fits in bins of (left) M^^o (right) cosO'^^y The shaded 
histogram shows the shapes expected from p^n^ Monte Carlo. 

parameters by ilcr); continuum rejection ±6.5% (by comparing the efficiency of the selection 
between data and MC for the -8+ D^n^ control sample); 7r° reconstruction efficiency 
±4.8% (by comparing the yields of 1] —>■ tt^tt^tt^ and — »• 77 between data and MC); 
tracking finding efficiency ±2.0% (from a study of partially reconstructed D* decays); pion 
identification efficiency ±0.8% (using the method described above). The contributions are 
summed in quadrature to obtain a total systematic error of 1}3;8%- 
We measure the branching fraction for —>■ p^vr^ to be 

B (5° ^ p^TT^) = (29.ll|;^(stat) ± 4.0(syst)) x 10"^ (1) 

A difference in the untagged decay rates to p^n^ and p^7r+ would indicate direct CP 
violation In order to measure this untagged asymmetry, we remove candidate events 
with ambiguous p charge assignment, i.e. events in the regions of the Dalitz plot where 
more than one of the combinations tt^tc^, tt~7c^, 7t~^tt~ are consistent with having originated 
from a p resonance. These are the regions where interference effects are strongest We 
fit the AE distributions for — * p+7r~ + B^ —>■ p+vr" (denoted a.s B —>■ p^tt~) and 
B^ —>■ p^7r+ + B^ p^7r+ {B p^7r+) candidates separately. The results are shown in 
Fig. El We find 36.7tlli B p+Tr" events, and 81.51^^:^ B p-n+ events. The charge 
asymmetry is calculated as 

^- NiB^ p+TT-) ± NiB ^ p-7r+) " U.^8_o.2i(statj_o.o5(systj, [2) 

where the systematic error is estimated by varying the PDFs used in the fit, allowing for 
asymmetry in the shape and normalization of the continuum component, and allowing for 
charge dependence in the efficiency. 

In order to obtain good sensitivity to the as yet unobserved decay 5° P^tt", additional 
discrimination against the continuum background is required. Therefore, a slightly different 
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FIG. 3: IS.E fit results for (left) B p^^n and (right) B ^ p ir^ . The components have the 
same meaning as Fig. ^ 

analysis procedure is followed. The majority of selection requirements are similar to those 
described above, however, neutral pion candidates are selected with a wider invariant mass 
window of 0.100 < My^/GeV/c^ < 0.165 to allow for the resolution of high momentum 
7r°s. One additional requirement is that possible contributions to the 7r"'"7r~7r° final state 
from b ^ c decays are explicitly vetoed. Another is that in order to reduce the event 
multiplicity before the best candidate selection, a smaller window in {AE, Mbc) is allowed. 
The selection region is -0.2 < AE/GeV < 0.4, Mbc > 5.23 GeV/c^. Whilst these changes 
do not significantly affect the p^n^ final state, they will be important in retaining a clean 
and unbiased 7r"'"7r~7r° Dalitz plot in the future. 

In order to select p°7r° from the three-body Tr+Tr^Tr" candidates, we require 0.50 < 
M^+^-/GeV/c2 < 1.10 and I cos^(^j,j| > 0.5. Contributions from B'^ — > p^n^ are explic- 
itly vetoed by rejecting candidates which fall into the invariant mass window of 0.50 < 
M^±^o/GeV/c^ < 1.10. This requirement also vetoes the region of the Dalitz plot where the 
interference between p resonances is strongest. 

The most notable difference from the selection requirements described above for the 
P^tt"^ final state is in the procedure to reject continuum background. Here, we make use 
of the additional discriminatory power provided by flavour tagging. In our published time- 
dependent analyses jl^], we define the variables q and r as being, respectively, the more 
likely flavour of the other B in the event {B^ {q = +1) or B^ {q = —1)), and a measure of 
the confidence that the q prediction is correct. As a corollary, events with a high value of 
\qr\ are well-tagged as either 5° or B^, and hence are unlikely to originate from continuum 
processes. Moreover, since the flavour tagging algorithm relies on particle identification 
information, it is unlikely that there is any strong correlation with any of the topological 
variables used above to separate signal from continuum. 

We use a large statistics sample of p^n^ MC and data from a continuum dominated 
sideband region (5.23 < Mbc/GeV/c^ < 5.26 and 0.2 < AE/GeV < 0.4), in order to 
simultaneously find the optimum selection requirements on \qr\ and Cg/ {Cg + Ch), as de- 
fined above. We find the sensitivity to p^n^ is maximized by requiring \qr\ > 0.74 and 
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FIG. 4: Effect of including \qr\ in the optimization procedure described in the text. The vari- 
ations in (top) S/\/S + B and (bottom) S/B with the \qr\ requirement are shown, after the 
Cs/ i^s + ^b) > 0.9 selection. The line indicates the selection requirement that is obtained from 
the optimization procedure {\qr\ > 0.74). Note that the statistical error on the entries becomes 
significant for values of \qr\ close to 1. For this reason, a smoothed histogram curve is also shown 
to guide the eye. 

Cs/ {^Cs + Ch) > 0.9, where the branching fraction of p°7r° is taken to be 1 x 10^^ as 

input to the optimization procedure. The effect of the \qr\ requirement is shown in Fig. |3 
Some modest improvement is seen in S/\/S + B, whilst S/B increases dramatically. 

As for p^Ti^, we obtain the signal yield by fitting the AE distribution after requiring 
events be in the Mbc signal region (Mbc > 5.269 GeV/c^), and cross-check the result by fitting 
the Mbc distribution after requiring events be in the AE signal region (—0.135 < AE/GeV < 
0.080). We model the signal using a Crystal Ball lineshape, with parameters determined from 
MC, and calibrated using control samples of B~ — > D^p~, with — > K~7t~^, p~ — >■ 7r~7r° 
and 5° — D*~^p~ , with D*^ — > D'^tt'^ and the same decays of the and p. In both cases 
the 7r° is required to have momentum in the CM frame greater than 1.8 GeV/c, in order to 
have similar dynamics to the p^Tc^ final state. We include additional components in the fit to 
accomodate continuum background (modelled by a Chebyshev polynomial with parameters 
obtained from a fit to the sideband region), p'^p^ (shape from smoothed histogram of MC 
events) and p'^Ti^ (shape from smoothed histogram of MC events). Contributions from b ^ c 
transitions are negligible in the fitted range of AE, whilst the tiny contributions possible 
from other rare B decays are accounted for in the systematic error. The only free parameters 
in the fit are again the signal and continuum normalizations; the p'^p'^ yield and polarization 
are fixed according to our recent measurements [l3|, while the p'^n^ yield is fixed based on 
theoretical expectations. 
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FIG. 5: Fits to (left) AE and (right) Mbc for p^i:^ candidates. In the AE fit, contributions from 
signal (dashed line) and p'^p^ (dot-dashed line) can be seen; the other contributions are very small. 
In the Mbc fit only contributions from signal (dashed line) and continuum are allowed. The sums 
of contributions are shown as solid lines. 



The fit result is shown in Fig. El The AE fit gives a signal yield of 6.6^2 g with a 
significance of 3.1a. The Mbc fit, in which only components for signal and continuum are 
included, gives a yield of 6.4114 with 3. 60" significance. 

In order to check that the signal candidates originate from — * p^n^ decays, and not 
from either non-resonant Tr+Tr^Tr" or crvr^, we relax the criteria on M^+^- and cos^^^j in turn 
and look at the distribution of the candidate events in both AE and Mbc signal regions 
in those variables. These distributions are shown in Fig. IHl Whilst the statistics are too 
small to make quantitative statements, there is no evidence for any contribution other than 
P^tt". We also consider possible contamination from p^yr^; from a large MC sample in 
which interference between resonances in not simulated, we expect < 0.1 events to pass the 
selection requirements. 

To obtain the branching fraction, we measure the efficiency using MC, and correct for 
the pion identification efficiency as above. The systematic error due to pion identification 
is ±3%. We also correct for possible differences between data and MC due to the \qr\ and 
Cs/ {Cs + Cb) selections; the statistical errors of the control samples {B^ D*^ p^ with 
D*+ D^T^+, K-TI+, p- 7r-7r° and B^ D^p~ with L)° ^ K-'k+, p- vr-vr" 

respectively) account for the largest contribution to the systematic error (±18%). We also 
calculate systematic errors due to: PDF shapes, by varying parameters by ±1ct (±3%); 7r° 
reconstruction efficiency, from the inclusive 77 study (±4%); track finding efficiency, from the 
partially reconstructed D* study (±2%). We repeat the fit after changing the normalization 
of the other B decay components according to the error in their branching fractions, and 
obtain systematic errors from the change in the result. For the unobserved mode p"'"7r° we 
vary the normalization by a factor of two. The total systematic eror due to rare B decays is 
±5%, we also verify that in the case that all rare B contributions are simulaneously increased 
to their maximum values, the statistical significance remains above 3a. The total systematic 
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error is ±20%, and we measure the branching fraction of —>■ p°7r° to be 

B (5° p°7r°) = (6.0l2|(stat) ± 1.2(syst)) x 10"^ (3) 

In order to test the robustness of this result, a number of cross-checks are performed. We 
vary the selection requirements on \qr\ and Cg/ i'^s + '^b)- In all cases, consistent central 
values for the branching fraction are obtained. We also select p^vr^ candidates after using 
this continuum rejection technique, and measure a branching fraction for p'^n^ which 

is consistent with that reported above. Furthermore, adopting the continuum rejection 
technique of our p^vr^ analysis and selecting p'^vr^ candidates also results in a consistent 
central value for the branching fraction of B^ P^tt*^, although the signal is insignificant 
above the large continuum background. 

In summary, we have measured the branching fraction 

B (5° ^ p^vr^) = (29.1+|;°(stat) ± 4.0(syst)) x 10"^ (4) 

in agreement with previous measurements. We also measure the untagged asymmetry to be 

^ = -0.38l°i?(stat)l°:°^(syst), (5) 

consistent with zero with the current statistical precision. In addition, we observe evidence, 
with 3.1(7 statistical significance, for B^ — > p'^vr^ with a branching fraction of 

B ^ p°7r°) = (6.0^|^(stat) ± 1.2(syst)) x 10"^ (6) 

This is the first evidence for B^ —>■ p^TC^, with a branching fraction higher than most 
predictions 0]. This may indicate that some contribution to the amplitude is larger than 
expected, which may complicate the extraction of 02 from time-dependent analysis of p^vr^. 
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